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NATIONAL  ADVISORY  COMMITTEE  FOR  AERONAUTICS 


A  TECHNICAL  NOTE  NO.  1543 


EFFECT  OF  CHORDWISE  LOCATION  OF  MAXIMUM  THICKNESS  ON  THE 
SUPERSONIC  WAVE  DRAG  OF  SWEPTBACK  WINGS 
By  Kenneth  Margolis 


SUMMARY 


On  the  "basis  of  the  linearized  theory  of  supersonic  flow,  equations 
are  derived  for  the  wave  drag  of  swept  "back  untapered  wings  at  zero  lift 
with  thin  double— wedge  sections  and  arbitrary  chordwise  location  of 
maximum  thickness.  Calculations  are  presented  for  a  representative 
supersonic  plan  form. 

The  optimum  location  of  maximum  thickness  for  untapered  supersonic 
wings  is  found  to  "be  at  50  percent  chord,  a  symmetrical  variation  in 
wing  wave-drag  coefficient  being  exhibited  about  this  minimum  value. 

The  drag  variation  is  slight  when  considerable  sweep  behind  the  Mach 
lines  is  present  and  the  variation  is  marked  at  Mach  numbers  where  the 
wing  approaches  or  is  swept  ahead  of  the  Mach  lines. 

It  is  found  tliat  for  tapered  plan  forms  the  sweep  of  the  line  of 
maximum  thickness  is  an  important  sweep  parameter  insofar  as  drag  due  to 
thickness  Is  concerned. 


INTRODUCTION 


A  recent  application  (reference  1)  of  the  linearized  theory  to  the 
calculation  of  supersonic  wave  drag  at  zero  lift  of  delta  wings  indicates 
that  minimum  drag  for  any  double-wedge  delta  wing  with  sufficient  sweep 
behind  the  Mach  lines  is  obtained  with  location  of  maximum  thickness  at 
10  to  20  percent  chord,  (in  reference  1  and  in  the  present  paper  the 
conventional  definition  of  "delta  plan  forms"  is  modified  to  include 
those  with  sweptback  trailing  edges.)  For  a  given  delta  plan  form,  a 
variation  in  maximum-thickness  location  necessarily  implies  a  variation 
in  the  sweep  of  the  line  of  maximum  thickness  and,  therefore,  no  general 
predictions  can  be  made  for  the  drag  effects  induced  by  each  parameter 
considered  separately.  In  fact,  there  Is  evidence  (see  reference  2)  to 
suggest  that  the  sweep  of  the  maximum-thickness  line  may  be  an  Important 
sweep  parameter  insofar  as  drag  due  to  thickness  is  concerned.  Hence, 
the  problem  of  Isolating  the  effects  of  chordwise  location  of  maximum 
thickness  immediately  presents  itself. 
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A 

This  problem  is  treated  in  the  present  paper  by  applying  the  method 
of  reference  3  to  derive  the  generalized  wave-drag  equations  of  swept back  . 

untapered  wings  at  zero  lift  having  thin  double-wedge  sections  with 
arbitrary  maximum-thickness  location.  The  wing  tips  are  chosen  parallel 
to  the  direction  of  flight  and  the  range  of  supersonic  Mach  number  for 
which  the  wing  is  swept  ahead  of  and  swept  behind  the  Mach  lines  is  con¬ 
sidered.  Typical  distributions  of  section  wave  drag  and  wing  wave-drag 
calculations  are  presented.  Comparisons  are  made  with  two-dimensional 
theory  and  with  results  obtained  for  tapered  wings. 


SYMBOLS 


x,y,z 

v 

p 

Ap 

<P 

M 


Cartesian  coordinates 
velocity  in  flight  direction 
density  of  air 
pressure  increment 
dynamic  pressure  (IpV2  ) 

disturbance-velocity  potential 
Mach  number 


dz/dx 

e 

c 

h 

t 

A 


m 


b 


slope  of  airfoil  surface 
half— wedge  angle 

chord  length,  measured  in  flight  direction 

location  of  maximum  thickness,  measured  from  leading  edge 

maximum  thickness  of  section 

angle  of  sweep,  degrees 

slope  of  leading  edge  (cot  A  ) 

span  of  wing 


Jt 


* 


D 


ving  area 
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A 


A 

X 


cd 


cd 


\ 

°®tip 

°D 


aspect  ratio  ^b^/S  y 

taper  ratio  (ratio  of  tip  chord  to  root  chord) 

section  wave-drag  coefficient  at  spanvise  station  y 
exclusive  of  tip  effect 

increment  in  section  wave— drag  coefficient  at  spanwise 
station  y  due  to  tip 

section  wave-drag  coefficient  at  spanwise  station  y 

(■Ho  +  °«tiP) 

wing  wave-drag  coefficient  exclusive  of  tip  effect 
increment  in  wing  wave-drag  coefficient  due  to  tip 
wing  wave-drag  coefficient  ( Cp  +  Cp^  j 


ANALYSIS 


V 


The  analysis  is  essentially  that  used  in  references  2  and  3.  For 
convenience,  a  brief  outline  of  the  basic  equations  is  presented. 

The  assumptions  of  small  disturbances  and  a  constancy  of  sonic 
velocity  throughout  the  fluid  lead  to  the  linearized  equation  for  the 
velocity  potential  cp 


(1) 


where  M  is  the  f>fech  number  of  the  flow  and  the  derivatives  are  taken 
with  respect  to  the  variables  x,  y,  and  z  of  the  Cartesian  coordinate 
system.  On  the  basis  of  this  linear  theory,  a  solution  for  a  uniform 
sweptback  line  of  sources  in  the  pressure  field  is  derived  in  reference  3. 
The  pressure  field  associated  with  this  solution  corresponds  to  that  over 
a  semi— infinite  sweptback  airfoil  of  wedge  section.  The  pressure  coef¬ 
ficient  Ap/q  at  a  spanwise  Btation  y  and  point  x  along  the  wedge 
for  m3  <  1  is 


Ap 

q 


—  =  —  tan  9 


m 


\I  1  —  m‘ 


2„2 


cosh 


— 1  x  —  m3 


Ply 


2 

y 

mx 


(2a) 


k 
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and.  for  m£  >  1 


^2  =  -  tan  0 

q  * 


m 


\/m^P^_—  1 


cos 


-1  X 


2 


P  y  -  mx| 


(2b) 


where  m  is  the  slope  of  the  leading  edge  of  the  wing,  0  Is  the  half- 


wedge  angle  (tan  0  ~  0  since  the  angle  is  small),  P  = 


M 


1  , 


and  the  origin  of  the  line  source  is  taken  at  ( 0,0).  In  the  region 
between  the  leading  edge  and  the  Mach  cone  (that  is,  ”■  X  <  y  <  mx), 

the  real  part  of  cos- ^  — ~~~  -mAZ  ig  constant  and  equal  to  it. 

Ply  -  mx| 

Equation  (2b)  then  reduces  to 


Ad 

q 


2  tan  0 


(2c) 


The  distribution  of  pressure  over  sweptback  wings  of  desired  plan 
form  and  profile  is  obtained  by  appropriate  superposition  of  wedge— type 
solutions.  In  order  to  satisfy  the  boundary  condition  over  the  surface 
of  an  untapered  wing  of  double— wedge  section,  semi— infinite  line  sources 
are  placed  at  the  leading  and  trailing  edge  of  the  wing  and  a  semi- 
infinite  line  sink  is  placed  along  the  line  of  maximum  thickness.  The 
strengths  of  these  lines  are  proportional  to  the  wedge  angles j  therefore, 
in  superimposing  solutions  of  the  type  given  In  equations  (2),  the 
appropriate  0  for  each  line  source  and  sink  must  be  found.  The  correct 
wedge  angles  necessary  to  obtain  a  double-wedge  section  of  length  c, 
maximum  thickness  t,  and  location  of  maximum  thickness  h  are  (see 
fig.  1) 
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^Leading 
line  source 


Trailing 
line  source 


Line  sink 


c 


2h 


c 

2(c  -  h) 


2h(c  -  h) 


\ 


►  (3) 

/ 


(It  is  convenient  to  express  the  angles  in  terms  of  the  section  thickness 

ratio  i  as  can  he  seen  by  reference  to  the  drag  integrals.)  At  the 

tip  where  the  wing  is  cut  off  in  the  flight  direction,  reversed  distri¬ 
butions  of  these  sinks  and  sources  are  placed  so  as  to  cancel  exactly  the 
effects  of  the  original  distribution  farther  spanwise  than  the  tip. 

Figure  1  shows  the  distributions  of  Bints  and  sources  for  an  untapered 
wing  of  double-wedge  section  and  Identifies  the  system  of  axes  and  symbols 
associated  with  the  derivation  of  the  drag  equations. 

The  disturbances  caused  by  the  elementary  line  sources  and  sinks  are 
limited  to  the  regions  enclosed  by  their  Mach  cones j  the  Mach  line  con¬ 
figuration  is  presented  in  figure  2.  For  purposes  of  simplification  the 
wings  considered  were  restricted  to  those  with  no  tip  effects  other  than 
the  effects  each  tip  exerts  on  its  own  half  of  the  wing. 

The  pressure  coefficients  obtained  from  superposition  of  solutions 
given  in  equations  (2)  are  converted  into  drag  coefficients  by  the 
following  relatione: 


2 

S 


cdc  ay 


Ap  dz 

q  dx 


dx 


dy 


(4) 


where  b  Is  the  wing  span;  S  is  the  wing  area;  dz/dx  is  the  slope 
of  the  airfoil  surface j  and  L.E.  and  T.E.  denote  leading  edge  and 
trailing  edge,  respectively. 


DERIVATION  OF  GENERALIZED  EQUATIONS 


The  drag  equations  are  derived  for  half  of  the  wing  since  the 
drag  is  distributed  symmetrically  over  both  halves.  The  induced  effects 
of  the  opposite  half— wing  are  represented  by  the  conjugate  terms  in  the 


where  A ,  B,  and  C,  which  refer  to  the  pressures  resulting  from  the  leading  line  sources,  line  sinks 
and  trailing  line  sources,  respectively,  are  as  follows: 
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where  D  and  E,  which  refer  to  the  pressures  resulting  from  the  leading 
line  sources  and  from  the  line  sinks,  respectively,  are  as  follows: 


D  =  cos  i 


2 

x  -  mp  y 
Ply  -  noc! 


+  cos 


O 

x  +  mg  y 

P|y  +  mx I 


cos 


_1  x  —  h  — 


mPfy 


p|y  -  m(x  -  h)| 


+  cos 


_1  x  -  h  + 


xnp2y 


P  |y  +  m(x  -  h)| 


It  should  he  noted  that  equations  (5)  give  the  drag  for  plen  forms  where 
the  tip  is  placed  farther  spanvise  than  the  points  of  intersection 
between  the  Mach  lines  and  the  leading  and  the  trailing  edges  (see  fig.  3). 
Deletion  of  certain  integrals  and  appropriate  changes  in  the  y— limits  of 
other  integrals  may  he  made  for  configurations  where  the  tip  is  placed 
nearer  the  root  chard.  Equations  (5)  are  evaluated  end  the  resulting 
section  wave-drag  and  wing  wave-drag  formulas  for  all  untapered  plan  forms 
are  nreeented  in  appendix  A. 

It  was  stated  previously  that  the  wings  considered  have  no  tip  effects 
other  than  those  each  tin  exerts  on  its  own  half  of  the  wing;  that  is,  the 
Mach  lines  from  one  tip  do  not  enclose  any  part  of  the  opposite  half— wing. 
This  condition  is  expressed  mathematically  as  fellows: 

For  c;8  <  1,  the  asnect  ratio 

A  >  j  (6a) 


and  for  mp>l,  the  aspect  ratio 


A  > 


2m 

mp'  +  1 


(6b) 


The  wavo-drag  contribution  of  the  tip  for  mP  <  1  is  (see  fig.  3) 
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TtS  /l  -  m2p2c 
8m(t/c)2  ^tip 


./I  ^2  pb/2 

*  v  1  —  mg  I  '  c  c 

hm(t/c)^  I  ad  ( 1+mP  )—mc  ^ip 

1+mB 


dy 


c2  rind 

Jmdfl+mP)— mh 
1+mP 


r£2* 

m  F  dx  dy 

Jd(l+mp)-yP 


Uh(c  —  h) 


rM 


f^+rcfO-mh 


+mp 


1+mp 


nXifflS 

m  Fdidy 

Jd(l+mp)-y3 


n  rad 

y+mc  1 

f  ffi— 

jmd(l+mP)-mh 

y+mh  F  te  ^ 

1+mP 

m 

i  md 

mdfl+mP  )-m(c  — h) 
1+mP 


y+mc 

m  G  dx  dy 
d(l+mp)+h-yp 


(T) 


where  F  and  G,  which  refer  to  the  pressures  resulting  from  the  leading 
line  sink  and  from  the  line  source,  respectively,  are  as  follows; 

F  =  cosh"1  K~ 

P|y  -  kx| 

G  =  cosh-1  *  r  ,.d  ~  A  t_2P2.(z -=JBgl 
P  y  -  m(x  -  h) 


For  mp>  1,  equation  (7)  is  still  valid  if  \l 1  —  m2P2 

m2P2  —  1  and  the  inverse  hyperbolic  function  cosh-1 

the  inverse  cosine  function  cos  1, 

Equation  (7)  is  evaluated  and  the  results  presented  in  appendix  B. 
The  total  wave-drag  coefficients  are  then  obtained  by  the  following 
relations: 


is  changed  to 
is  changed  to 
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Cd  =  +  C\ip 

> 

Op  =  Op  + 
v  v*>  tip 


(8) 


It  is  found  that  0 


’D. 


tip 


is  identically  equal  to  zero  and,  hence. 


Cp  =  Cj^  for  wings  satisfying  the  aspect— ratio  limitations  imposed  in 
equations  (6) . 


When  the  expression  for  Cp  is  differentiated  with  respect  to  h 
in  order  to  find  an  optimum  maximum-thickness  location,  it  is  found  that 


(9) 


and  that  the  drag  function  is  minimized  at  this  value  of  h  independently 
of  the  other  parameters.  It  is  also  seen  freon  the  drag  equations  that 
the  distributions  are  symmetrical  about  this  minimum  point. 


RESULTS  AND  DISCUSSION 


For  calculation  purposes,  an  untapered  plan  form  of  aspect  ratio  2 
and  sweepback  of  60°  has  been  used.  Equation  (9)  indicates  that  the 
trends  obtained  for  this  representative  plan  form  are  equally  valid  for 
all  untapered  wings  with  double-wedge  profiles.  The  results  may  be 
assumed  to  have  qualitative  application  to  curved  profiles  without  cusps. 

Section  wave  drag.—  Spanwise  distributions  of  section  wave-drag 
coefficient  are  shown  in  figures  b  and  5  far  Mach  numbers  of  l.kl^  and 
respectively.  At  a  Mach  number  of  l.hlk  the  wing  is  swept  behind  the 
Mach  lines,  and  at  a  Mach  number  of  3  the  wing  is  swept  ahead  of  the 
Mach  lines.  In  each  figure  maximum-thickness  locations  are  varied  from 
20  percent  chord  to  80  percent  chord.  Variations  in  maximum-thickness 
location  need  actually  only  be  considered  up  to  50  percent  chord;  the 
results  for  any  arbitrary  percent;  chord  k  and  (100  —  k)  are  equal 
because  of  symmetry  considerations.  (See  .drag  equations  in  appendix  A.) 
When  the  wing  is  swept  behind  the  Mach  lines  (fig.  *0,  the  centroid  of 
the  drag  forces  on  a  wing  panel  moves  noticeably  inboard  with  forward 
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or  rearward  shift  of  the  maximum  thickness  from  the  midpoint  position. 
When  the  wing  is  swept  ahead  of  the  Mach  lines  (fig.  5)>  little  change 
is  noted  in  the  centroid. 

Wing  wave  drag.—  Variation  of  wing  wave-drag  coefficient  with  Mach 
number  is  shown  in  figure  6.  As  was  mentioned  previously  in  the 
mathematical  discussion  of  the  drag  equations,  the  symmetrical  section 
is  seen  to  be  the  optimum  profile  for  an  untapered  plan  form  insofar  as 
minimum  drag  is  concerned.  The  variation  of  wing  vave-drag  coefficient 
with  maximum-thickness  location  does  not  appear  to  be  too  significant 
when  the  wing  is  swept  far  behind  the  Mach  lines;  however,  as  the 
"critical"  condition  is  approached  (that  is,  m£  =  l)  the  variation 
becomes  more  noticeable.  There  is  seen  to  be  a  marked  drag  variation 
with  maximum-thickness  location  in  this  region  and  for  Mach  numbers  where 
the  wing  is  swept  ahead  of  the  Mach  lines. 

Figure  7  is*  in  effect,  a  cross  plot  of  figure  6  and  presents  vari¬ 
ation  of  wing  wave-drag  coefficient  with  maximum-thickness  location  for 
Mach  numbers  of  l.kl4  and  3.  The  previous  remarks  are  clearly  illustrated 
in  this  figure  by  the  flat  curve  at  the  lower  Mach  number  and  by  the  curve 
with  a  well-defined  minimum  at  the  higher  Mach  number. 

Comparison  with  two-dimens i  nnal  theory.  —  If  m  approaches  infinity 
(that  is,  a  rectangular  plan  fom),  the  drag  equation  reduces  to 


CD  =  1  cg 

(t/c)g  ^Mg  —  1  k(c  ~  h) 


(10) 


which  is  exactly  the  result  obtained  by  Ackeret  for  two-dimensional 
flow  (see  reference  1).  This  result  is  expected  since  the  plan  forms 
considered  have  zero  increment  in  wing  wave  drag  due  to  the  wing  tip. 

The  two-dimensional  (Ackeret  solution)  and  three-dimensional  (swept- 
back  wing)  results,  therefore,  exhibit  the  same  qualitative  drag  effects 
due  to  variation  in  maximum-thickness  location  —  that  is,  a  drag  vari¬ 
ation  symmetrical  about  a  minimum  value  at  50  percent  chord.  Of  course 
the  quantitative  results  are  quite  different,  the  three-dimensional 
value  being  lower  when  considerable  sweep  behind  the  Nkch  lines  is  present 
and  the  two-dimensional  value  being  less  for  other  conditions.  Figure  8 
presents  the  variation  of  two-dimensional  drag  with  Mach  number  for 
different  maximum-thickness  locations.  Comparison  with  the  three- 
dimensional  result  for  the  wing  of  aspect  ratio  2  and  sweepback  of  60° 
is  also  indicated  in  this  figure. 


Comparison  with  tapered  wings.—  The  results  presented  in  reference  1 
for  a  given  delta  plan  form  (X  =  0)  indicate  an  optimum  maximum¬ 
thickness  location  at  10  to  20  percent  chord  when  the  wing  is  swept 
sufficiently  behind  the  Mach  lines.  The  results  obtained  in  the  present 
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investigation  for  untapered  wings  (X  =  1)  indicate  an  optimum  maximum- 
thickness  location  at  50  percent  chord.  Hence,  the  optimum  location  for 
any  arbitrary  and  conventional  tapered  wing  (0  S  H  l)  with  sufficient 
sweep  behind  the  Mach  lines  approaches  the  10-percent  to  20-percent  value 
for  small  taper  ratios  (large  taper)  and  the  50— percent  value  for  large 
taper  ratios  (small  taper). 

For  a  given  delta  plan  form,  the  sweep  of  the  maximum-thickness  line 
varies  with  the  location  of  maximum  thickness  and  thus  no  prediction  can 
be  made  for  the  drag  effects  induced  by  each  parameter  considered  separ- 
rately.  For  an  untapered  wing,  however,  the  effects  of  maximum-thickness 
location  are  isolated  since  the  sweep  of  the  maximum-thickness  line 
remains  constant.  Inasmuch  as  the  present  results  indicate  an  adverse 
drag  effect  due  to  shifting  the  maximum-thickness  location  forward  of  the 
50-percent  location,  the  optimum  location  of  maximum  thickness  for  the 
delta  wings  of  reference  1  of  between  10  and  20  percent  chord  must  result 
from  the  increased  sweep  of  the  line  of  maximum  thickness  at  the  forward 
location.  The  sweep  of  the  line  of  maximum  thickness  thus  appears  to  be 
an  important  sweep  parameter  for  tapered  wings  insofar  as  drag  due  to 
thickness  is  concerned.  This  inference  is  further  supported  in  reference  2 
in  which  the  section  wave-drag  coefficient  at  the  root  of  tapered  wings  is 
found  to  be  a  function  of  the  Mach  number  and  the  sweep  of  the  maximum¬ 
thickness  line  and  is  found  to  be  independent  of  leading-edge  and  trailing- 
edge  sweep. 


CONCLUSIONS 


The  following  conclusions  refer  specifically  to  nonlifting  wings 
with  double-wedge  profiles  but  may  be  assumed  to  have  application  to 
curved  supersonic  profiles  without  cusps: 

1.  The  optimum  location  of  maximum  thickness  for  untapered  supersonic 
wings  is  at  50  percent  chord,  a  symmetrical  variation  in  wing  wave-drag 
coefficient  being  exhibited  about  this  minimum  value  similar  to  that  found 
in  two-dimensional  supersonic  flow. 

2.  The  variation  of  wing  wave— drag  coefficient  with  maximum-thickness 
location  is  Blight  when  the  wing  is  swept  far  behind  the  Mach  lines  and 

is  marked  at  Mach  numbers  where  the  wing  approaches  or  is  swept  ahead  of 
the  Mach  lines. 
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3.  The  sweep  of  the  line  of  maximum  thickness  is  an  important 
sweep  parameter  for  tapered  wings  insofar  as  drag  due  to  thickness  is 
concerned. 


Langley  Memorial  Aeronautical  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Langley  Field,  Va.,  October  31#  19^7 
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APPENDIX  A 

EVALUATION  OF  EQUATIONS  (5)  FOR  SECTION  WAVE  DRAG  AND  WING 
WAVE  DRAG  EXCLUSIVE  OF. TIP  EFFECTS j 
Section  Drag  for  m3  <  1 

\ll  -  m2p2 


IAVE  DRAG  A] 

5^  <  0.5^ 


4m(t/c)2 


c  -  A 


for 


o<  y<  j 


mh 


m3 


=  A  +  B 


=  A  +  B  +  C 


=  A  +  B  +  C  +  D 


mh 
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Figure  1.-  Symbols,  wedge  angles,  and  distributions  of  sinks  and  sources  for 

an  untapered  wing. 


Figure  2.-  Mach  line  configurations  for  untapered  plan  forms. 


Figure  3.-  Information  pertinent  to  integration  limits  in  equations  (5)  and  (7) 
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Figure  6.-  Variation  of  wing  wave -drag  coefficient  with  Mach  number  for  various 
maximum -thickness  locations.  Aspect  ratio,  2;  taper  ratio,  1;  sweepback 
angle,  60°. 


Figure  7.-  Variation  of  wing  wave -drag  coefficient  with  maximum 


thickness  location.  Taper  ratio,  1;  aspect  ratio,  2;  sweepback 
angle,  60°. 


Figure  8.-  Two-dimensional  variation  of  wing  wave -drag  coefficient 
with  Mach  number  for  various  maximum -thickness  locations  and 
comparison  with  three-dimensional  result  for  untapered  wing. 
Aspect  ratio,  2;  sweepback  angle,  60°. 


